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Abstract 

We study the critical behavior for inhomogeneous versions of the Curie-Weiss model, where the 
coupling constant J^- (d) for the edge ij on the complete graph is given by J^j (/3) = Wk)- 

We call the product form of these couplings the rank-1 inhomogeneous Curie-Weiss model. This model 
also arises (with inverse temperature /3 replaced by sinh(/3)) from the annealed Ising model on the 
generalized random graph. We assume that the vertex weights (wi)ie[Af] ^-re regular, in the sense that 
their empirical distribution converges and the second moment converges as well. 

We identify the critical temperatures and exponents for these models, as well as a non-classical 
limit theorem for the total spin at the critical point. These depend sensitively on the number of finite 
moments of the weight distribution. When the fourth moment of the weight distribution converges, then 
the critical behavior is the same as on the (homogeneous) Curie-Weiss model, so that the inhomogeneity 
is weak. When the fourth moment of the weights converges to infinity, and the weights satisfy an 
asymptotic power law with exponent r with t S (3,5), then the critical exponents depend sensitively 
on r. In addition, at criticality, the total spin Spf satisfies that /7 VC-i)/R- 2) converges in law to 
some limiting random variable whose distribution we explicitly characterize. 


1 Introduction 

Universality is a key concept in the theory of phase transitions, with application to a large variety of 
physical systems. Informally, universality means that in the thermodynamic limit different systems show 
common properties close to criticality. The theory based on the renormalization group suggests that 
systems fall into universality classes, defined by the values of their critical exponents describing the nature 
of the singularities of measurable thermodynamic quantities at the critical point. 

In the presence of heterogeneities, e.g. spin systems on random graphs used to model interaction on 
a network m da US 121] it is not clear a-priori to what extent universality applies. From the point of 
view of the structure of the network, emerging properties of real networks have been identified in several 
empirical studies in different contexts - social, information, technological and biological networks. Many 
of them are scale free, with a degree sequence obeying power-law distribution, and small world, with short 
graph-distance among vertices. As a consequence power-law random graphs, i.e., graph sequences where 
the fraction of nodes that have k neighbors is proportional to k~'^ for some r > 1, are often used as 
mathematical models for real-world networks. In this paper we investigate universality for spin system on 
power-law random graphs displaying phase transitions. 
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The issue of universality is related to the network functionality. Indeed the occurrence of a thermody¬ 
namic phase transition is associated to a change in macroscopic properties of the networks, for instance 
the possibility to reach consensus in a social network can be related to the occurrence of a spontaneous 
magnetization. Thus the investigation of differerrt universality classes for spin systems on random graphs 
is a relevant question with immediate practical relevance for the network functionality. 

Due to the random environment, when considering the Ising model on the random graphs used to 
model real networks, a distinction is required between different averaging procedures. Two settings are 
often studied in the literature: the quenched measure (graph realizations are studied one-by-one so that 
they produce a random Boltzman-Gibbs measure) and the annealed measure (all graph realizations are 
considered at once and they give rise to a deterministic Boltzman-Gibbs measure). See |18l [7] for an 
extended discussion of the two settings. 

In the paper m the quenched critical exponents have been rigorously analyzed for a large class of 
random graph models. More precisely in im it is proved that the critical exponent d (describing the 
behavior of the magnetization at the critical temperature as the external field vanishes), the exponent 
f3 (describing the behavior of the spontaneous magnetization as the temperature increases to the critical 
temperature) and the exponent 7 (describing the divergence of the susceptibility as the temperature 
decreases to the critical temperature) take the same values as the mean-field Gurie-Weiss model whenever 
the degree distribution has a finite fourth moment. This includes for instance the case of the Erdos-Renyi 
random graph. For power law random graphs, it is proved that for r > 5 the model is in the mean-field 
universality class, whereas the critical exponents are different from the mean-field values for r € (3,5). 

In this paper we provide the analysis of the critical behavior but in the annealed setting. Our results are 
fully compatible with the universality conjecture. The annealed critical temperature is different (actually 
higher) than the quenched critical temperature, but the set of annealed critical exponents that can be 
rigorously studied are the same as the quenched critical exponents. In the annealed setting our results are 
stronger since we are able to show that 7 ^ = 7 . Here 7 ' describes the divergence of the susceptibility as 
the temperature approaches to the critical temperature from below, and in the quenched setting we were 
able to show only that 7 ^ > 7 . 

A main difference between the quenched and annealed case is that while the analysis of the quenched 
measure could be done in great generality, the study of the annealed case is much harder. Indeed the results 
of HH are valid for all graph sequences that are locally like a homogeneous random tree [210 n Em | 22 ] 
and uniformly sparse. For the annealed setting it is not enough to control the behavior of the model on the 
typical graph realizations (namely rooted random trees). For the annealed measure one needs to study 
exponential functionals of the graphs, i.e., questions on large deviations of sparse random graphs that 
are largely unsolved. Thus we specialize our analysis of the annealed critical exponents to a particular 
class of random graphs models. This is given by the Generalized Random Graph models, also called 
inhomogeneous random graphs of rank -1 in the literature (see [20111] for a non-rigorous study). 

By exploiting the factorization of the Gibbs measure and the edges independence we reduce the study 
of the annealed measure for the Ising model on the Generalized Random Graph to the analysis of an 
inhomogeneous Curie-Weiss model. As we shall see, for this model we are able to also study the properties 
at criticality. On a sequence of temperatures approaching the critical value, we prove the scaling limit for 
the properly renormalized total spin. As a result, our findings extend the analysis of the scaling limit of the 
standard Curie-Weiss model [laESlEl] and provides new asymptotic laws for the (properly renormalized) 
total spin. 
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2 Model definitions and results 


2.1 Inhomogeneous Curie-Weiss model 

We start by defining the inhomogeneous Curie-Weiss model. This is a generalization of the classical Curie- 
Weiss model in which the strength of the ferromagnetic interaction between spins is not spatially uniform. 
As the standard Curie-Weiss model, it is defined on the complete graph with vertex set [N] := {1,..., A^}. 
See Table 1 at the end of the paper for a summary of the important notation used in this paper. 

Definition 2.1 (Inhomogeneous Curie-Weiss model). Let a = {cdieiAr] £ variables. 

The inhomogeneous Curie-Weiss model, denoted by CWjv(T), is defined by the Boltzmann-Gibbs measure 

gi^jv(o-) 

Tn{(t) = —— ( 2 . 1 ) 

Zjn 

where the Hamiltonian is 

= \ + B Yi ( 2 - 2 ) 

i,je[N] is [AT] 

and Zjv is the normalizing partition function. Here j3 is the inverse temperature, B is the external magnetic 
field and J = are the spin couplings. 

In the above, the interactions Jij{l3) might be arbitrary functions of the inverse temperature (in 
particular no translation invariance is required), provided that the thermodynamic limit is well-defined, 
i.e., the following limit defining the pressure exists and is finite, 

(j){l3,B) := lim ^ log Zj.,(13, B) . (2.3) 

N^oo 

In the following we will restrict to the ferromagnetic version of the model, i.e., we will assume Jij{(3)>0. 
Since the coupling constants J = {Jij{f3)}ij^[N] ^^e positive and possibly different for different edges, we 
speak of an inhomogeneous Curie-Weiss model. We next state our hypotheses on the coupling variables. 
Each vertex i G [N] receives a weight Wi, We will take J = {Jiji/3)}ij£[N] such that 

where ^ w'fc- (2.4) 

^ k£[n] 

In the case where Wi = 1, our model reduces to the (homogeneous) Curie-Weiss model. We will call the 
coupling constants in (12.41) the rank-1 inhomogeneous Curie-Weiss model. In Section 12.31 we describe the 
assumptions that we make on the weight sequence w = (rci)jg[ 7 v]- 

In |18| it is shown that the rank-1 inhomogeneous Curie-Weiss model arises in the study of the annealed 
Ising model with network of interactions given by the rank-1 inhomogeneous random graph, also called 
the generalized random graph, which we describe next. 

2.2 Generalized random graph 

In the generalized random graph mm, each vertex z G [A^] receives a weight Wi. Given the weights, edges 
are present independently, but the occupation probabilities for different edges are not identical, rather 
they are moderated by the weights of the vertices. We assume that the weights w = (rci)jg[ 7 v] are strictly 
positive (there is no loss of generality in supposing this, since the vertices with zero weight will be isolated 
and can be removed from the network). 
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Definition 2.2 (Generalized random graph). Denote by Iij the Bernoulli indicator that the edge between 
vertex i and vertex j is present and by pij = P {lij = 1) the edge probability, where different edges are 
present independently. Then, the generalized random graph with vertex set [iV], denoted by GRGjv('ie), is 
defined by 

WiWi 

pb = 1 )—^- 

tjv + WiWj 

where t'jv = total weight of all vertices. 

We have now defined two classes of models that depend on vertex weights w = (rCi)jg[Ar]. We next 
state the assumptions on these weights. 


2.3 Assumptions on the vertex weights 

We study sequences of inhomogeneous Curie-Weiss models and generalized random graphs as N —)• oo. For 
this, we need to assume that the vertex weight sequences w = are sufficiently nicely behaved. Let 

Upf G [N] denote a uniformly chosen vertex in GRGjv('U^) and W)v = wjJm weight. Then, the following 
condition defines the asymptotic weight W and set the convergence properties of (VFiv)iv>i to W: 

Condition 2.3 (Weight regularity). There exists a random variable W such that, as N ^ oo, 

(i) A w, 

fnj E[W2] = ± ZielN] ^ E[W2] < oo, 
where —?• denotes convergence in distribution. Further, we assume that E[LF] > 0. 

Note that, by uniform integrability, Condition I2.3r ii) implies that also E[VF)v] = E[VF] < 

oo. 

Condition 12.31 implies that the sequence (GRGjv('m))Ar>i is a uniformly sparse tree-like graph with 
strongly finite mean and with asymptotic degree D distributed as a mixed Poisson random variable. 


P(D = /c) = E 



vpfc- 

TT 


( 2 . 6 ) 


see e.g., m Chapter 6]. 

Our results depend sensitively on whether the fourth moment of W is finite. When this is not the 
case, then we will assume a power-law bound on the tail of the asymptotic weight: 

Condition 2.4 (Tail of LF). The random variable W satisfies either of the following: 

(i) E[W4] < oo, 

(ii) W obeys a power law with exponent r G (3,5], i.e., there exist constants > 0 and tco > 1 

such that 

< P(W > w) < \/w > Wo. (2.7) 


To prove the results on the scaling limit at criticality we will strengthen 

Condition 2.5 (Tail of Wpp and deterministic sequences). The sequence 
either of the following: 

fi) E[W4] = A wf ^ E[W4] < oo, 


our assumptions as follows: 
of weights (rCi)ig[Ar] satisfies 
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(a) it coincides with the deterministic sequence 


Wi = C^{-r) ’ 

for some constant > 0 and r G (3, 5). 

We remark that the above deterministic sequence is A^-dependent (we do not make this dependence 
explicit) and its limit W satisfies (|2.7I) since Wi = [1 — F]~^{i/N), where F{x) = 1 — for 

w > Cu)- In the next section, we explain what the annealed measure of the Ising model on GRGjv('U^) is. 


2.4 Annealed Ising Model 

We hrst dehne the annealed Ising model in general on finite graphs with N vertices, then we specialize 
to GUGpf^w). We denote by Gjv = (V)v,I?jv) a random graph with vertex set Vjv = [A] and edge set 
Eff GVff X Vff. We denote by the law of the graphs with N vertices. 


Definition 2.6 (Annealed Ising measure). For spin variables a = (cri, taking values on the space 

of spin configurations fijv = { — 1,1}'^ the annealed Ising measure is defined by 


PNia) = 


(exp B Eie[iV] ) 


Q^{Zg^ (/3,R)) 


where 

fd ^ (Tiaj + R ^ (Tj 

is the partition function. 


Zgn iP, P)= exp 


(2.9) 


( 2 . 10 ) 


With abuse of notation in the following we use the same symbol to denote both a measure and the 
corresponding expectation. 


Definition 2.7 (Annealed thermodynamic quantities). For a given A G N we introduce the following 


thermodynamics quantities at finite volume: 


(i) The annealed pressure; 

M(d,B) = llog(Q^(Z^(/3,i?))). 

(2.11) 

(ii) The annealed magnetization; 


(2.12) 

where the total spin is defined 

as 



Sn = . 

(2.13) 


ie[N] 


(Hi) The annealed susceptibility; 

XNif3,B) = Am^(/3,R). 

(2.14) 
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2.5 Annealed Ising Model on GRG 


We now specialize the previous definitions to the annealed Ising Model on the Generalized Random Graph. 
By assuming the probability pij of each edge in is that given in (12.Sp . we can compute explicitly the 
average of the partition function (12.101) . Indeed, recalling that /jj is the indicator of the edge between 
vertex i and j, we can write 


Qn (Zn ^ exp 1^/3 ^/ijCTiCTj + R ^ cjj 

o-eOjv i<j ie[iv] 

and, by using the independence of the variables /jj, we compute [18] that 

where C'(/3) > 0 is a constant and the positive couplings Jij{l3) are defined as 

1 , f e^pij + {1 - Pij) 




^Pij + (1 - Pij) 


(2.15) 


(2.16) 


(2.17) 


The r.h.s. of (|2.16p can be seen as the partition function of an inhomogeneous Gurie-Weiss model with 
couplings J given by (I2.17p . Thus, the annealed Ising model on the GRGjv('*a’) is equivalent to such 
GW^(J), i.e., the two measures coincide point-wise on the sample space. Our proof (see eq. (I4.63p l shows 
that the Jij{/3) in (I2.17p are close to the form in (j2.4p with /3 replaced by sinh(/3), so that the study of 
the annealed generalized random graph reduces to the rank-1 ICW model. Preliminarily to the statement 
of our main results we recall the model solution given in m- By symmetry, we always take B > 0. We 
denote by (3c the annealed critical inverse temperature defined as 


/3c := inf{/3 > 0 : M(/3,0+) > 0}, (2.18) 

where the spontaneous magnetization is given by 

M(/3,0+)= lim lim Mjv(/3,B) . (2.19) 

B^0+ N^oo 

Theorem 2.8 (Thermodynamic limit for annealed Ising on GRGjv('if) and for rank-1 CWjv(T) |18|). Let 
(Gjv)jv>i be a sequence of GIiGr^{w) graphs satisfying Condition \2.3[ Then, 

(i) For all 0 < (3 < oo and for all i? € R, the annealed pressure exists in the thermodynamic limit 
A —>■ oo and is given by 

V^(/3,R) := lim V'^(/3,R). (2.20) 


(ii) The magnetization per vertex exists in the limit A —)• oo and is given by 

M{P,B) := lim M^{f3,B). 

N—>-oo 


( 2 . 21 ) 


The limit value M equals: M{(3, B) = -^ifi/d, B) for B>0, whereas M = 0 in the region 0 < /3 < j3c, 
R = 0. More explicitly, when B > 0 or B = 0'^ and (3 > /3c 


M{/3,B) = E 


tanh 


/sinh (/3) 
E[iy] 


Wz* + B 


where z* = z*{/3,B) is the unique positive solution of the fixed point equation 

z = E tanh 


I sinh {(3) 


E[IT] 


Wz + B 


I sinh (/3) 


E[IT] 


W 


( 2 . 22 ) 


(2.23) 
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(in) The annealed critical inverse temperature is given by 


where 


Pc = asinh {l/n) , 

_ E[VK2j 

^ “ E[VF] ■ 


(iv) The thermodynamic limit of susceptibility exists and is given by 


X{P,B) 


lim Xn{P,B) 
N^oo 


92 

dB^ 


H(3,B). 


(2.24) 

(2.25) 


(2.26) 


(v) For the rank-1 inhomogeneous Curie-Weiss model CW{^{J), (i)-(iv) hold with P replaced with asmh{P). 

Theorem 12.81 shows that a phase transition exists for the annealed Ising model on the generalized 
random graph and the rank-1 inhomogeneous Curie-Weiss model. For the rank-1 inhomogeneous Curie- 
Weiss model in the special case where Wi = 1, Theorem 12.81 reproves the classical result for the Curie-Weiss 
model. When the weights are inhomogeneous, the critical value is instead given by Pc = 1/v. 

Let us compare the annealed critical value in (12.241) to that in the quenched setting as derived in [17| . 
There, it is proved that the quenched critical value Pc'^ equals P'p^ = atanh(l/j^) > asinh(l/z^) = Pc- Thus, 
the annealed critical value is smaller due to a collaboration of the Ising model and the graph properties. 

In this paper, we analyze the block spin scaling limits at Pc and we study the universality class of 
the model. For this, we define the annealed critical exponents analogous to the random quenched critical 
exponents as in HU: 

Definition 2.9 (Annealed critical exponents). The annealed critical exponents /3, 7 , 7 ' are defined by: 


MiP,0+)^iP-Pcf, 

for P \ Pc] 

(2.27) 

M{Pc,B)^B^I\ 

for B \0] 

(2.28) 

x(/3,0+)x(/3,-/3)-^. 

for P / Pc] 

(2.29) 

x(/3,0+)x(/3-/3e)-^'. 

for P \ Pc, 

(2.30) 


where we write f{x) g{x) if the ratio f{x)/g{x) is bounded away from 0 and infinity for the specified 
limit. 

We remark that, as is customary in the literature, we use the same letter for the inverse temperature P 
and one of the magnetization critical exponent f3. In this paper they are distinguished by the use of the 
plain, respectively bold, character. 


2.6 Main results 

We start by proving that the annealed critical exponents for the magnetization and the susceptibility take 
the values conjectured in |20| . 

Theorem 2.10 (Annealed critical exponents). Let (Giv)Ar>i be a sequence o/GRGjv(rr) graphs fulfilling 
Conditions \2.‘J[ and \2.4\ Then, the annealed critical exponents defined in Definition \2.y\ using Pc given in 
(|2.24p exist and satisfy 
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re (3,5) 

E[W^] < 00 

/3 

l/(r - 3) 

1/2 

S 

r-2 

3 

7 = 7' 

1 

1 


For the boundary case r = 5 there are the following logarithmic corrections for (3 = lj2 and S = 3; 

The same results hold for the rank-1 inhomogeneous Curie-Weiss model CW^(J), the critical exponents 
being now defined using j3c = Ijv. 


Remark 2.11 (Comparison to the Curie-Weiss model). For the rank-1 inhomogeneous Curie-Weiss model, 
we see that the inhomogeneity does not change the critical behavior when the fourth moment of the weight 
distribution remains finite, but it does when the fourth moment of the weight distribution increases to 
infinity. In the latter case, we call the inhomogeneity relevant. 

Remark 2.12 (Comparison to the quenched case). In fll^ . the first two and fourth authors of this paper 
have shown that the same critical exponents hold for the quenched setting of the Ising model on power-law 
random graphs, such as GRGjv('U^), under the assumptions in Conditions and \2.4\ In m, however, 
we only managed to prove a one-sided bound on 7 '. Thus, our results show that for GRGjv(tf) both the 
annealed and quenched Ising model have the same critical exponents, but a different critical value. This is 
a strong example of universality. 

Remark 2.13 (Extension of 7 = 1). The result 7 = 1 holds under more general conditions, i.e., E[1E^] < 
00 . See Theorem \3. 61 below. 

From the previous theorem we can also derive the joint scaling of the magnetization as {/3,B) \ (/3c, 0): 
Corollary 2.14 (Joint scaling in B and (/? — fdc))- For t 5, 

MiffB) = e{{P-fIcf + BF^), (2.32) 


where f{/3,B) = Q{g{l3,B)) means that there exist constants ci,Ci > 0 such that cig{f3,B) < f{ffB) < 
Cig{ll, B) for all B € (0, e) and /? G (/3c, jdc + £) with e small enough. For r = 5, 


M{/3,B) = 0 


l3-(3c y/2 / B y/3N 

logl/(^-^c)^ Vlog(l/R)/ / 


(2.33) 


Our second main result concerns the scaling limit at criticality. The next theorem provides the correct 
scaling and the limit distribution of Sjv at criticality (for a heuristic derivation of the scaling, see the 
discussion in Section [2.711 . For GRGiv('ir’), we define the inverse temperature sequence 


/3c, AT = asinh(l/i/jv), 


(2.34) 


where 

_eM 

" E[W^]’ 

so that /3c,AT /3c for N ^ 00 . For rank-1 CW^iJ), we replace /3 by asinh(/3), so that /3c ,n 

main result is the following: 


(2.35) 


I/vn- Our 
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Theorem 2.15 (Non-classical limit theorem at criticality). Let (Gjv)Ar>i be a sequence of GRGjv(k;) 
graphs satisfying Conditions \2.3\ and Condition 1^.51 and let S have the respective values stated in Theo- 
rem \2.1(K Then, there exists a random variable X such that 


S, 


V 


Ar<5/(5+i) 


X, 


as N ^ oo, 


(2.36) 


where the convergence is w.r.t. the measure Pm at inverse temperature fic^N = asinh(l/i/jv) and external 
field B = 0. The random variable X has a density proportional to exp(—/(x)) with 


fix) = < 


1 nw^] 4 

12 E\WT^ 


E,: 


1 f t-2 


when < oo, 

when r G (3, 5). 


(2.37) 


The same result holds for the rank-1 inhomogeneous Curie-Weiss model at its critical value fic^N = 


We will see that in both the case where the fourth moment is finite as well as when it is infinite, 


lim 44 = C, 


a:—>oo X 


l+<5 


with 


C = 


1 ¥.[Wf 
12 E[WR 

r—1 


(2.38) 


(2.39) 


when E[W^] < oo, 

— log cosh (ly when r G (3,5). 

This result extends the non-classical limit theorem for the Gurie-Weiss model to the annealed GRGjv(if) 
and the rank-1 CWjv(T). 


2.7 Discussion 

Random weights. Instead of choosing the weights w deterministically, one can also choose the weights 
i.i.d. according to some random variable W, with E[IT^] < oo. In this case. Condition 12.31 holds a.s. by 
the laws of large numbers. Hence, if Qm denotes the average over all graphs drawn according to the GRG 
conditioned on the weights, then our results also hold a.s. When in the annealing also the average over the 
weights is taken, then the model becomes unphysical, because the pressure becomes infinite as is proved 
in dH]. 

Critical exponents. Theorem 12. lOl imnlies that the annealed exponents are the same as in the quenched 
case. Indeed, by (12.6p . the condition E(IT'^) < oo is equivalent to E(iG^) < oo, where K is the forward 
degree of the branching process describing the local structure of GRGjv(rc). Thus the conditions in 
Theorem 12.101 defining the universality classes are the same as those in Theorem 2.8 in |11| . 


Scaling limit of block spin variable. In |18| . it is proved that the classical central limit theorem for 
the total spin Sm holds in the one-phase region of the annealed Ising model i.e.. 


Sm — PNiSn) V ^ 

y/N 


W(0,x), w.r.t. PjV) as N ^ oo. 


(2.40) 


In |17| we prove the analogous result in the quenched setting. More precisely, we prove (|2.40p for the 
quenched measure in the quenched uniqueness regime for all random graphs that are locally tree-like. A 
prominent example is the GRG]v('U^) as studied here. 
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At criticality, i.e. for (/3, B) = (Pc, 0), the limit in (12.401) is no longer true. A scaling different from ^/N 
has to be used to obtain the scaling limit, and also this limit is not a normal random variable. In |15[ I16| . 
Ellis and Newman prove that for the standard Curie-Weiss model 


S, 


V 


iV3/4 




as N ■ 


oo, 


(2.41) 


where X is a random variable with density proportional to exp{ —We extend this result to the 
rank-1 inhomogeneous Curie-Weiss model, and thus to the annealed Ising model. We prove that the 
scaling with is also correct when E[1E^] < oo, but different for r G (3,5). Furthermore we show that 
when = oo, different asymptotic distributions arise in the scaling limit. We characterize them for 

the weight deterministic sequence (|2.8p in which the weights follows a precise power-law. Such a sequence 
is rather generic in the sense that it produces an asymptotic weight that is also power-law distributed. 
The analysis shows that the fluctuations of the total spin decrease as the exponent r becomes smaller and 
the distribution seen in the scaling limit has tails proportional to e“^*^ . 


Heuristic for the scaling limit. To obtain a guess for the correct scaling, we can use the standard 
scaling relation between S and ij as in M- On a box in the d-dimensional lattice with side lengths n, 
[n]*^ C the exponent r] satisfies 

(2.42) 


where Pn^'^ is the expectation w.r.t. the Ising measure on this box and Sn is the sum of all spins inside the 
box, where it should be noted that there are sites in the box. Hence, to compare this with our setting, 
we take N = and, with an abuse of notation, let Sn = •S'jv- If there is an exponent A such that 
converges in distribution to a non-trivial limit, then it must also hold that Pjv(('S'n/-Al'^)^) = PN^Sn/n^'^^) 
converges. Hence 5^ ~ so that d -|- 2 — r; = 2dA. The standard scaling relation 2 — rj = d|^ |14) 
now suggests that we should choose 



(2.43) 


We prove that this is indeed the correct scaling and we also show that the tail of the density behaves like 
exp{—as is conjectured on (see m Section V. 8 ]). 


Near-critical scaling window. Theorem 12.151 is proved along the critical sequence /3c ,at approaching 
the critical inverse temperature /3c in the limit N ^ oo. A different scaling limit might be obtained by 
working with a sequence near the critical one, the so-called near-critical window, i.e., /3c ,at -|- Ajv with 
—>• 0 at an appropriate rate. As is argued in Section [4.51 it turns out that for the annealed Ising model 
the width Ajv of the scaling window is and the scaling limit differs by a quadratic term that 

appears in in the function f{x) describing the density of in (I2.37p . 


At criticality. As a consequence of the previous discussion, we also infer that if one works at critical 
inverse temperature /3c, the scaling limit that will be seen to depend on the speed at which i^jv approaches 
n. Indeed, from (I2.24j) and (I2.34p . one has /3c — I3c,n = 0{v — v^). For a natural example given by 
the deterministic sequence in Condition 12.51 (ii) one has that when r > 5 then u — = o(l/iV^/^) 

and thus the limiting distribution does not change; on the contrary when r G (3, 5] then ly — = 

-|- o(l)) for some 7 ^ 0 , and thus the distribution changes since we are shifted in the 
near-critical window. See again Section 14.51 for more details. 


10 














Organisation of this paper. In Section [3l we start by deriving the annealed critical exponents in 
Theorem 12.101 In Section 21 we prove our non-classical limit theorems at criticality in Theorem 12.151 We 
will prove our results only for the annealed GRGjv('IA’), since the proofs for the rank -1 inhomogeneous 
Gurie-Weiss models are either identical, or simpler. 


3 Annealed critical exponents: proof of Theorem 12.10 

We follow a strategy similar to that in m, although the proof in our case is a bit easier since the annealed 
magnetization is expressed in terms of the deterministic fixed point z* in (I2.23p . whereas in the quenched 
setting the magnetization is expressed in terms of a fixed point of a distributional recursion. The proof of 
Theorem 12.101 is split into Theorems 13.51 dealing with the exponents (3 and 6 (Section 13.ip . Theorem 13.61 
for the exponent 7 and Theorem 13.71 for the exponent 7 ' fSection 13.2p . Some lemmas and propositions 
containing preliminary results are also stated and proved in Section 13.11 

Our analysis of the critical behavior crucially builds on the fixed point equation (|2.23p . We apply 
truncation arguments together with monotonicity (see the proof of Proposition 13.31 for a prototypical 
example). We rely on Taylor expansion properties for the fixed point z* in (|2.23p as is customary for the 
Ising model. By truncation we mean that we decompose the range on integration of various expectations 
with respect to the limiting distribution W according to the size of the fixed point z* and using asymptotics 
for truncated moments of W. 


3.1 Magnetization: critical exponents f3 and 6 

We start by showing that the phase transition is continuous. 

Lemma 3.1 (Gontinuous phase transition). Let be a sequence with /?£ and non-increasing, 

Pi > Pc CLnd Bi > 0, and Pi \ Pc and Bi\^Q as i ^ 00 . Then, the solution of (12.231) satisfies 


lim z*{Pi, Be) = 0. (3.1) 

£—>00 

In particular, 

\\m. z*{Pc,B) = D, and lim z*{P,D'^) = Q. (3-2) 

Proof. The existence of the limit (j3.1jl is a consequence of the monotonicity of z*{P,B) and the fact that 
z*{P,B) > 0 for B > 0. Suppose that limi^ao z*{Pe, Be) = c > 0. Then, it follows from (12.231) and 
dominated convergence that 


c = lim z*{Pi,Bi) = E 

£^00 


tanh 


/ / sinh (Pc) 

IE[IP] 



sinh (Pc) 
E[IT] 


W 


< csinh(Pc)r' = c, 


(3.3) 


where we used that tanh(x) < x for x > 0 and Pc = asinh(l/i/). This contradiction proves the lemma. □ 


We next show that z* has the same scaling as we want to prove for M(P,B) by proving the upper and 
lower bounds in Propositions 13.31 and 13.41 below. These then allow us to obtain the theorem. But first we 
state some properties for truncated moments of W in the following lemma: 


Lemma 3.2 (Truncated moments of W). Assume that W obeys a power law for some r > 1, see item 
(a) in Condition\2.4\ Then there exist constants Ca,T,Ca,T > 0 such that, as £ ^ 00 , 


Ca,rr-(-l) 
Cr-1,T log^ 


< E [W“]l{^<,}] < 


n pa-(T-l) 

Cr-l,r log^ 


when a > T — 1, 
when a = T — 1. 


(3.4) 
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and, when a < t — 1, 


(3.5) 

□ 


Proof. The proof is similar to that of HH Lemma 3.4]. 

In the following we write Cj, Ci, i > 1 for constants that only depend on fd and on moments of W and 
satisfy 

0 < liminfcj(/?) < limsupci(/3) < oo, (3.6) 

and the same holds for Ci. The constants Ci appear in upper bounds and Cj in lower bounds. Furthermore, 
we write e*, i > 1 for error functions that depend on (5, B and on moments of W, and satisfy 

limsupei(/3,i?) < oo and lim ei(/3c, i?) = 0. (3-7) 

B\o 

Here, the subscript i is just a label for constants and error functions. 

Further, we introduce the following notation that will be used extensively in the following: 

- Vw- 

Proposition 3.3 (Upper bound on z*). Let j3 > fdc and B > 0. Then, there exists a Ci > 0 such that 

z* < -^£[1^] sinh(/3)H + sm\i{l3)vz* — Ciz*^, (3.9) 

where S takes the values as stated in Theorem \2.1(K For t = 5, 

z* < ■\/E[IU] sinh(/3)i? + sinh(/3)z/2;* — Ciz*^ log {1/z*). (3.10) 

Proof. We frequently use that tanh(H) < B. A Taylor expansion around x = 0 gives that, for some 

C G (0,x), 

tanh(x + B) = tanh(i?) + (1 — tanh^(H))x — tanh(H)(l — tanh^(H))x^ 

— -^(1 — tanh^(C + B))x^ + tanh(C + i?)(l — tanh^(C + B))x^ 

O 

1 Q 4 o 

< B + X — -x^ + - tanh(x + B)x^, (3-11) 

O O 

where we also used that tanh(x) <1. If we now assume that x + B < atanh|, then 

tanh(x + B) < B + X — -x^. (3.12) 

6 

We apply this result to (|2.23p where x = a{l3)Wz*, which we force to be at most atanh| by introducing 
an indicator function as follows: 


z* < E [{B + a{P)Wz*) a(/3) IF] 


+ E 


{tanh {aif3)Wz* + B) - {B + a{P)Wz*)} a(/3) 


(3.13) 


since tanh(i? + x) < B + x. Hence, using (I3.12p . 


2- < VE|M']sinh(^))B + sinh(/J)^2- - -a(/ 3 )*E 




(3.14) 
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t-5 


For E[VF^] < oo, this is indeed of the form (|3.9I) and we are done. If r G (3,5), then it follows from 
Lemma 13.21 that 

4 ] f ct(/S) 

^ l{a(/3)r^.*+S<atanhi}J > ^4,^ (^^anhi - i?) ' 
which proves the proposition for r G (3,5). The proof for r = 5 is similar and we omit it. □ 


E 


(3.15) 


We now proceed with the lower bound: 

Proposition 3.4 (Lower bound on z*). Let 13 > Pc and B > 0. Then, there exists a ci > 0 such that 

z* > VE[IT] sinh(/3)S + sinh(/3)i/ z* - - Bei, (3.16) 

where S takes the values as stated in Theorem \2.1(A For r = 5, 

(3.17) 


z* > y^E[IF] sinh(/3)i? + sui\i{P)vz* — ciz*'^ log {l/z*) — Bei. 
Proof. As in (|3.11jl we can bound 

tanh(x + B) > B + X — — B{B + Bx + x^), 

O 


(3.18) 


where we have used that B — B^ < tanh(i?) < B. For E[IF^] < oo, we can immediately use this to obtain 


where 


and 


2 ;* > -^E[IF] sinh(/?)i? + smh.{P)uz* — Ciz*^ — Bei, 
1 E[IF^] 


Cl = -sinm(/3) 


E[W] 


2 ’ 


ei=B^E [W] sinh (/3) + B smh{P)iyz* + 


(3.19) 

(3.20) 

(3.21) 


All terms in ei indeed converge to 0 in the appropriate limit, because of Lemma [34 
For r G (3,5), we rewrite 2 ;* as 


z* = \/E[IF] sinh(/3)i? + smh.{P)vz* 

+ E [{tanh {a{P)Wz* + B) - {B + a{P)Wz*)] a{P) W (liw<i/.*} + hw>i/z*})] ■ (3-22) 

The case where W < 1/z* can be treated as above. This gives 


E 


tanh {a{P)Wz* + B) - {B + a{P)Wz*) \a{P) Wl{w<i/z*} 


> -^sinh2(/3) 


E[WH{w<i/z>}\ 


z*^-Be2, 


where 

62 = .Bq;(/3)E[IF]1{i^<i/^*}] + Bsmli{P) 
By Lemma 13.21 


E[IF]2 


z* + a{PrE[W^l[W<i/z*}]z*^- 


E[IF] 

E[W%w<i/.*}\ < C^,rz*^-^ 


(3.23) 

(3.24) 

(3.25) 
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so that indeed (13.231) is bounded from below by 


with 


- C2Z 


—2 


C 2 = ^ sinh^(/3) 


Be2. 

(3.26) 

C4,r 

E[IU]2‘ 

(3.27) 


Using Lemma’s 13.II and 13.21 one can also show that all terms in 62 indeed converge to 0 in the appropriate 
limit. 

It remains to bound the term where W > Xjz*. For this we use that tanh(x + i?) > 0: 


E 


where 

By Lemma [3.21 


tanh {a{l3)Wz* + B) - {B + a{(3)Wz*) |a(/3) Wl{w>i/z*} 

-Eippi-^ -Be3, 

63 = a(/ 3 )lE[LF]l{vi/>i/^*}]- 


(3.28) 

(3.29) 


nw%w>i/z>}] < C2,rZ*^-\ (3.30) 

again giving the right scaling. As a consequence (13.281) is bounded from below by — Be^ with 

C 3 = ^ sinh(/3) 


Similarly, 

63 < a(/S) 

which indeed converges to 0. We conclude that (I3.16P holds with ci = C2 + C3 and ei = 62 + 63. □ 


CI 2 ,. 

E[IU] 2 ’ 

(3.31) 

*r-2 

(3.32) 


The upper and lower bounds on z* in the previous two propositions allow us to prove that the critical 
exponents take the values stated in Theorem 12.101 

Theorem 3.5 (Values of /3 and S). The critical exponents [3 and 6 equal the values as stated in Theorem 
when E[LF^] < 00 and r G (3, 5). Furthermore, for r = 5 (I2.3ip holds. 

Proof. Proof for exponent (3. We start by giving upper bounds on the magnetization. From (I2.22p it 
follows that 

M(/3, B) = E [tanh {a{P)Wz* + B)] < B + ^E [LF] sinh (/3)z*. (3.33) 

We first analyze (3 and hence take the limit B 0 for (3 > /3c. This gives 

M(/3,0+) < VE[LF]sinh(/3)zo, (3.34) 

where we write Zq = lims^o -z*- Since M(l3, 0"*“) > 0 by the definition of /3c, the same must be true for Zq. 
We will deal first with the cases r € (3,5) and E[VF^] < 00 . Taking the limit B \ 0 in (13.Op and dividing 
by Zq, we get for r / 5 

CiZq^~^ < sinh(/3)zz — 1, (3.35) 

and hence, observing that P = 1/{S — 1), 

Zq < (sinh(/3)i^ — 1)^ . (3.36) 
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From a Taylor expansion of sinh(/3) around Pc = asinh(IE[VF]/E[VF^]) it follows that 

sinh(/3)z^ — 1 < cosh(/3)i/(/3 — /3c). (3.37) 

Hence, 

M{P, 0+) < v'E [VF] sinh (cosh(/3)jy)^ (/3 - Pcf, (3.38) 

so that it is easy to see that 

r M{P,0+) ^ 

limsupy-- < oo. (3.39) 

/3\/3c (P - Per 

The lower bound can be obtained in a similar fashion. Starting from tanh x > x — and taking the limit 
H 0 for /3 > /3c in (|2.22p . we obtain 

M(/3,0^) > -\/E \W] sinh {P)zq — sinh(/3)j^Zo^. 


Again, starting from the lower bound (I3.16p . taking H 0 and dividing by Zq 

Zq > (sinh(/3)i/ — 1)^ , 

and, by a Taylor expansion around /3c, 

sinh(,d)z/ - 1 = cosh(/3c)z/(/3 - /3c) + 0((/3 - /3c)^). 

Using (j3.3Hp . fl3.4ip and (I3.42P in (I3.40h we obtain: 

M(/3,0+) > Vie [VF] sinh(/3)c7^ [cosh(/3c)z/(/3 - /3c) + 0((^ - ^c)")] ^ 


(3.40) 

(3.41) 

(3.42) 


- sinh(/3)i/C'i [cosh(/3c)z^(/3 - /3c) + 0((/3 - /3c)^)]^^ , 


2 m 2 P 


which shows that also 


0 < liminf 


A7(/3,0+) 


(3.43) 


(3.44) 


/3\/3c (/3-/3c)P’ 

concluding the proof for the exponent /3 in the cases r G (3,5) and E[VF^] < oo . In the case r = 5 we 
can prove the upper bound for M(/3,0'*') in a similar fashion, i.e., taking the limit H 0 for /3 > /3c of 
(I3.10p and dividing by Zq. This yields to 


^*2 ^ sinh(/3)i/ - 1 ^ cosh(/3)z/(/3 - /3c) ^ ~_(/3 - /3c) 


Ci\og{l/z*Q) C'ilog(l/V) 


log(lVS)’ 


(3.45) 


where (|3.37p has been used in order to obtain the second inequality and cosh(/3) has been bounded in a 
right neighborhood of /3c to obtain the third inequality. Since x i—>■ l/log(l/x) is increasing in (0,1) and 
Zq < C{P — Pc) 2 for some C > 00 form (I3.45P we obtain: 


2*2 < 
Zq ^ 


C{P-Pc 


Cilog(l/[C(/3-/3c)V2])- 

The previous inequality together with (13.341) . proves the upper bound 

M(/3,0+) 


lim sup — , 

( p-Pc ^ ^ 

llog(l/(/3-/3c)) 


< oo. 


(3.46) 


(3.47) 


^The proof of Zg < U//? — /3c) ^ can be obtained by rewriting (13.4511 as —Zq^ log < k{/3 — Pc), for some k > 0. Since 
w < —wlogta for w < 1/e, we conclude that for /3 — /3c > 0 small enough, the previous inequality gives Zq'^ < k{P — Pc)- 
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The lower bound can be obtained in the same way. Indeed, from (13.171) in the limit B ^ 0, we obtain, for 
some positive constants C and C 


2 ^ sinh(^)z^ - 1 ^ p, (13- Pc) 


7 > 

Zq ^ 


> c 


> c- 


{P-Pc 


Ci\og{l/zl) \og{l/zl) \og{l/{P - Pc)y 


(3.48) 


where, once more, we have used that x e-)■ l/log(l/x) is increasing in (0,1) and the bound Zq > C{P — 
for some C > 0 and any 0 < £ < 2^ The previous inequality plugged in (I3.40p gives 


( 13-13, 


1/2 


> 0 , 


(3.49) 


\^log(l/(/3-^c)), 

concluding the proof for r = 5. □ 

Proof for exponent 5. We continue with the analysis for <5. Setting /3 = /3c in (13.9p . we obtain 

z*{Pc,B)<{Ci^W^])-"^By\ (3.50) 

Using this inequality in (13.331) with P = Pc, we obtain 


which proves that 


M{Pc, B) <B+ " ^i/<^ 

(yEiU^)i+i 


V M{Pc,B) ^ 

lim sup- —n — < oo 

B\0 


since <5 > 1. Inequality (I3.16p with P = Pc gives 




(3.51) 

(3.52) 

(3.53) 


This estimate, along with (I3.50p . will be used in the lower bound of the magnetization at /3 = /3c obtained 
by tanhx > x — 


M{Pc,B) >B + 


E[IU] 

v/e[U^ 


:(1 - 2B)z*{Pc, B) - z*{Pc,Bf - By 


giving, for B > 0 small. 


1/5 


(3.54) 


M(/3c, B)> B+ (1 - 2B)c~^^^ ( ^ ^ - eiiPc, B) 1 B^/^ 

- (Ci7E[IU])"2/^52/^ - B^. 

Recalling that lim^^o Ci(/3c, 73) = 0 and 5 > 1, the previous bound gives 

limigf > 0, 


(3.55) 


(3.56) 


which concludes the proof for S in the cases r G (3,5) and E[IU^] < 00 . The analysis for r = 5 can be 
performed in a similar way as for (3. □ 


Proof of Corollary \2.14\ The proof can be simply adapted as in [111 Corollary 2.9]. □ 

^The proof of tlie inequality Zq > C{P — for 0 < e < 2 can be obtained starting from the rightmost inequality 

of (13.481) combined with the fact that log 1/a; < AeX~^ for all x £ (0,1) and any e > 0. 
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3.2 Susceptibility: critical exponents 7 and 7 ' 

We now analyze the susceptibility and compute the critical exponents 7 and 7 '. We start by computing 
the former under more general conditions than those of Theorem 12.101 


Theorem 3.6 (Value of 7 ). For E[1T^] < 00 , 

^Im^ x(/3, 0+)(/3e - /?) = tanh(^,), 

so that 7 = 1 . 

Proof. From Theorem 12.81 it follows that in the one-phase region, i.e., for /3 < /3c or 7 ^ 0, 


(3.57) 


xiP,B) = ^M{/3,B)=E 


8z* \ 

+ a{P)W— j (1 - tanh2 (a(j9)lVz* + B)) 


We can also compute the derivative of z* by taking the derivative of (I2.23P : 


dz* 

'dB 


= E 


,dz* 


a{P)W + aiPfW^— ) (1 - tanh^ {a{l3)Wz* + B)) 


(3.58) 


(3.59) 


If we take the limit B 0 for /3 < /3c, then the tanh^(-) term vanishes, since by definition of /3c it holds 
that Zq = lim^x^o -2^* = 0- Hence, if we write 


then (I3.59P simplifies to 


Solving for ^ gives 


^ = lim -^z*{P,B), 
dB B\odB 


r)7* r)?* 

= VE[W] sinh(/3) + sinh(/3)i.^. 


dzl _ v^E[IT] sinh(/3) 


dB 1 — sinh(/3)z^ 

Also taking the limit H 0 in (|3.58p and using the above gives 

1 — smh.[p)i' 

From a Taylor expansion around /3c, we get that 

sinh(/3c) - cosh(/3c)(/3c - /3) < sinh(/3) < sinh(/3c) - cosh(/3)(/3c - /3), 

so that 

E[IT]2 sinh(/3) 

1 + < x(/3,0+) < 1 + 


(3.60) 

(3.61) 

(3.62) 

(3.63) 

(3.64) 


E[VF]2sinh(/3) 


E[IT2] cosh(/3)(/3c-/3)^ 


(3.65) 

□ 


E[IT2] cosh(/3c)(/3c -/3) 
since sinh(/3c)z^ = 1. Hence, (I3.57P follows. 

We now analyze 'j': 

Theorem 3.7 (Value of 7 '). For W satisfying Condition \2.4\ with E[IT^] < 00 or with r G (3,5), 

7 ' = 1. (3.66) 
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Proof. We split the proof into the two cases that cover the hypotheses of the theorem. 

(a) Proof under the assumption E[VP^] < oo. We are now in the regime where /3 > /3c, so that Zq > 0. 
We start from (|3.59p . take the limit i3 \ 0 and linearize the hyperbolic tangent. In order to control this 
approximation, we define g{x) = — tanh^(x) and remark that on the basis of our assumption on W, we 

have that E[(iy^ V l)g{W)] < oo. It will be useful also to factorize g{x) = x‘^k{x) with k{x) = 0(1) as 
X —>■ 0, so that we also have E[W®A:(iy)] < oo. This gives 


r)7* 

^ =E 
dB 


r)7* \ 

(/3)W + a{(3)V-^j (1 - tanh^ (a(/3)WzS)) 


= VlE[IP] sinh(/3) - eo + ^ (sinh(/3)z^ - a(/3)^ ¥.[W^]zf) + 


dB 


+ E 


dz. 


-^a{l5fW^ + aim ) 9 WW4) 


where 


Solving (|3.67|) for ^ gives 


/ sinh(/3) 
V E[IT] 


3/2 

E [W^] zf. 


dzl _ sinh(/3) - eo - E [a(/3)Wff {aim4)\ 

dB 1 - sinh(/3)jz + a{(5Y ^ {aim^)]' 


(3.67) 

(3.68) 

(3.69) 


To analyze (I3.69P we use the lower and upper bounds in Propositions 13.41 and 13.31 Taking the limit B 7^0 
in ()3.16l) with d = 3, ci given in (|3.20p and dividing by Zq, we obtain 


> 3 *=1'‘T 


(sinh(/3)zz — 1). 


E[IT^] sinh^(/3) 

Taking the same limit i? 0 in (|3.14l) and dividing by Zq we obtain also 

E[IT]2 1 


< 6 


lE[^^l{„(/3)M/^*<atanhi}] sinh2(^) 


(sinh(/3)z/ — 1). 


(3.70) 


(3.71) 


By Taylor expansion, 

sinh(/3)z^ - 1 = cosh(/3c)zz(/3 - /3c) + 0((/3 - /3c)^), (3.72) 

we conclude, from (j3.70p . (j3.7ip . and the fact that ^\W^^^a{^p^\!Yz*<SLta.nh^\ 1E[IT^] as /3 —>■ /3c, that = 

0(/3 — /3c). Using this, we can now evaluate the terms in numerator and denominator of (I3.69h as /3 ^ /3c. 
The first term in the numerator has a non vanishing finite limit, while eo = 0(/3 — /3c). The third term (ig¬ 
noring the irrelevant multiplicative factor a(/3) ) is 'E\Wg {aif3 )WZq)] = a(/3)^2;Q^E \yV^k iaifi)WZ q)\ = 
0((/3 - /3c)). Indeed, since fc(x) < 


a(/3)So*^ E [W^k [aim4)] < = 0(/3 - /3c). (3.73) 


Let us now consider the denominator and define 


Z?(/3) := 1 - sinh(/3)zz + a(/3)^ ElW^- (3-74) 

By (I3T(1 . (13TB and fiTTB . 

2 cosh(/3c)i/(/3 - /3c) + 0((/3 - ^ < Zl(/3) < (a(/3) - 1) cosh(/3c)i/(/3 - /3c) + 0((/3 - m, (3-75) 
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where a(/3) is a function that converges to 6 as /? —>■ /3c- Thus, from the previous display we obtain 
T)(/3) = 0(/3 — /3c). The fourth term in the denominator of (|3.69p . again discarding an irrelevant factor 
and arguing as before, is E \W‘^g Zq)\ = a(/3)^2;Q^E \W^k {a{l3 )WZq)\ = — f5cf‘). Therefore, 

summarizing our hndings, 

r)z* 

^=0((/3-ft)-‘). (3.76) 

From (I3.58p . the upper bound follows using (|3.76l) : 


x(/3,0) < E 




< 1 + Vsinh(/3)E[bF]0((/3 - (3^)-^). 


(3.77) 


Similarly, for the lower bound we use that 1 — tanh^(x) > 1 — for every x, we obtain 


x(/3,0) > E 


.dzX 




= 1 +E 


dB 

a(/3)M^|| -E[a(/3)2fyV] -E 




= 1 + Vsinh(/3)E[bF]0((/3 - /3c)-^) - sinh(^)z/0(/3 - /3c) - a{pfK[W^]0{l), (3.78) 


again starting from (|3.58l] . using (|3.76p and = 0(/3 — /3c). From (|3.77l) and (|3.78p we obtain 

0 < liminfx(^,0+)(/3-/3c) < limsupx(/3,0+)(/3 -/3c) < oo, (3.79) 

/3\/3c /3\/3c 

proving the theorem in the case that E[VF^] < oo. 

(b) Proof for W satisfying Condition 12.41 (ii). Now we generalize the previous proof in order to 
encompass also the case of those W whose distribution function F{w) = 1 — P(VF > w) satisfies Condition 
I2.4r iih We start by defining 

hg^B,z* (w) = tanh {awz* + B) aw — a^w‘^z*, (3.80) 

where the dependence of a on /3 has been dropped, and rewriting (12.231) as 

z* = E + a22*E[w2]. (3.81) 


Using integration by parts, 

/* + CxD /• + CXD 

^[hp,B,z*iW)] = / hp^B,z*{w)dF{w) = - hi3^B,z*iw)d{l - F{w)) 

Jo Jo 

= - lim [hi3 B,z*{w){l - F{w))] + hp^B,z*{0){l - F{0)) 

UI^+OO 


+ 


r+oo 

/ h'i^,B,z*iw){l - F{w))dw. 

Jo 


The boundary terms in the previous display vanish and therefore 

f+OO 


E 


r+oo 

[h0,B,z*{W)]= / h'p^B,z*{w){l- F{w))dw. 

Jo 


(3.82) 


(3.83) 


Taking into account that the power law of Condition I2.4f iil holds for w > wq, we write the previous 
integral as 

E [hp^B,z* (W)] = G(/3, B,z*) + J(/3, B,z*), (3.84) 
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where 


rWQ /• + OD 

G{/3,B,z*):= h'p^B,z*iw){'i- - Fiw))dw, J{P,B,z*):= - F{w))dw. (3.85) 

Jo j WQ 


Therefore, (|3.8ip can be rewritten as 

2 * = G(/3, B, z*) + J(;0, B,z*)+ 


(3.86) 


Now we take the limit i? 0 in the previous equation. Recalling that Zq := lim^^o z* > 0) since the 
following limits exist: 


Vim G{l3,B,z*) = G{l3,Zo), lim J{/3,B,z*) = J{j3,ZQ) 

-D\0 i? \0 


(3.87) 


by bounded convergence, then we arrive to 

= G(/3, zS) + J{P, z*o) + (3,88) 

In the next step we bound J{/3,Zq). From the dehnition of J{/3,B,z*) in (|3.85p . and Condition I2.4r iil. 


/• + CXD /• + CXD 

Cw / h'p^Q.,*{w)w~^'^~^^dw<J{j3,B,z*)<Cyv / 

JwQ J WQ 


(3.89) 


Applying the change of variable y = az*w leads to 


p-\-oo r-\-oo 

/ = [taiih{y + B) — ytanh‘^{y + B) — y] y~^'^~^'^dy. (3.90) 

JwQ JaWQZ* 


Therefore, denoting 


r+oo 


I{/3,B,z*) := / [tanh(y+ R)-ytanh^(y+ R)-y] y ^Uy, 


• aWQZ^ 


we can rewrite ^3.891) as follows: 

c^a^-h*"-‘^I{^,B,z*) < J{/3,B,z*) < Gy,a^-^z*"-‘^I{P,B,z*). 
Since, again by bounded convergence. 


^+oo 

hm/(/3,R,z*) = / [tanh(y) - ytanh2(y) - y] =:/(/3,Zo), 

JcxWQZq 


B\0 

we obtain from (I3.89P that 

c^a^-^zl^-^I{P,zl) < Ji/3,z*o) < G^cF-^zl^-^I{l5,zl). 
On the other hand, since tanh(y) — y tanh^(y) — y < 0 for y > 0, we also have 

r+cxD 


(3.91) 


(3.92) 


(3.93) 


(3.94) 


/ -hOO 

[ytanh2(y) + y -tanh(y)]y"(^"^)dy < -/(/3,Zo) 

r+oo 

< / [ytanh^(y)+y-tanh(y)]y"('^"^)(iy =: iF(T). 

Jii 


(3.95) 
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Therefore, from (|3.88p . (|3.94p and (13.951) 


Zq > G{l3,Zo) - cwa^~^Zo^~^K{T)+a^ZoK[W'^] 


(3.96) 


and 


z*o < G{P, z*q) - Gwa^~^z*o^-‘^k{T) + a^z*QE[W'^] 


(3.97) 


The next step is to control the behaviour of G{I3,Zq) as (3 ^ 13c- We start by showing that G(/3,2 ;q) is 
O(zq^) as 13 ^ (3c- From the definition of G{(3, Zq), 


rwQ 

G{I3,Zq)= / [—— tanh^(azot(;)a^Zotc + atanh(azot(;)] (1 — T(t(;))dt(;. 

Jo 

Since the function between the square brackets is negative for y > 0 and decreasing, we have 

j-Wo 

0 >G{(3,Zq) > [—o^wqZq — tanh^(Q;ri;o-2o) + atanh(amo^;o)] / (1 — F{w))d 

Jo 


(3.98) 


> [—a^WQZQ — a'^WQZQ ianlr{aw^ZQ) + a tanh(at(;o5;o)] = --a^W^ZQ + 0{zq) 

O 


(3.99) 


where the last equality is obtained by Taylor expansion. 

Thus, the previous inequality implies that G(/3,2 ;q) = O(zq^). Again, from (I3.96|) and (|3.97p dividing 
by zl, 

1 - a‘^E[W^] > z*^-^ (G(/3, z^)zf-^ - c^a-^-^Kir )), (3.100) 

and 

1 - a2E[lF2] < ^*r-3 zl)zf-^ - Gwa^-^k{T )). (3.101) 

Since G{I3-,Zq)z(^~'^ = 0{z^~'^) and r G (3,5), the previous inequalities together with (I3.72p imply that 
z*o^~^ = 0{I3 - 13c) as ^ \ /3c- 

Next, we consider the derivative of Zq. Again, taking the limit B '\0 for (3 > f3c of (13.591) we obtain 


dz^ aE[W]-aE[Wtanh2(aWz())] 

^ ~ 1 - a2E[lF2] + a2E[lT2 ia-nh^{aWzl)]' 


(3.102) 


Since the numerator has a finite positive limit as (3 (3c (in particular, the second term is vanishing), we 

will focus on the denominator 

D2{P) := 1 - a‘^E[W^] + a^E[W^tanh.^{aWzo)]- (3.103) 

We start by decomposing the average 

E[W^tanh^(aWzQ)] = E[W^ tanh^/aWzQ/t^w^wo}]^anh.^{aWzo)l{w>wo}], (3.104) 
and analyze the two terms separately. The first one can be bounded as follows 

0 < E[1T^ tanh2(aVF2;Q)]lny<^Q}] < (3.105) 

showing that 

E[WHanh\aWz*o)l{w<wo}] = 0{zf) = Oi{(3-/3c)^), (3.106) 
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with the exponent satisfying 2/(r — 3) > 1 since r G (3,5). The second term can be treated with the 
integration by parts formula 


E[lT^tanh^(Q;fT 2 :Q)]lr^>^ i] = — lim tanh^(ati; 2 :Q)(l — ^(tc))] 

UI^ + OO 


+ 


r+oo Q 

ujq tanh^(arco2;o)(l — ^(uJo)) + / [rc^ tanh^(Q!r(;zo)](l — F(t(;))dt(;. (3.107) 

J OW 


Since r > 3, from Condition 12.41 we conclude that the limit in the previous display vanishes. It is also 
simple to see that 

wl taiih^{awoZQ){l - F{wo)) = 0{zf) = 0{{/3 - /3c)~). 

From (j3.1fl3p and using (I3.104p . (I3.106p . (|3.107p . (13.108^ . we can write 


Zl2(/3) = ^2(/3) + 0((/3-/3c)^), 


(3.108) 

(3.109) 


with 


D2iP) := (l-a2E[VF2]) + a- 


r+oo 

/ 


dw 


[up'ia-n}c?{awZ q)]{ 1 — F{w))dw. 


(3.110) 


The second term in the r.h.s. of (j3.109p is 0{{j3 — /3c)®) with s > 1, therefore we can forget it since the 
first term of D2{(3) is 0(/3 — /3c). Now we focus on the second term of D2{/3). 

By using (|2.7I) and applying the change of variable y = az*w, we can bound the integral in the last 
display as 

r-+oo Q 

[w'^ ianlP{awz’P)\{l — F{w))dw < CwcF~^ M{t), (3.111) 


/ 

J WO 


dw 


where 


f + CxD 


M{t) := / [2y tanh^(7/) + 2?/^ tanh(7/)(l — tanh^(y))] y 


(3.112) 


and the bound in (I3.11ip is obtained thanks to the positivity of the integrand. The convergence of the 
integral is ensured by the fact that this function is 0{y^~'^) close to y = 0 with 1 > 4 — r > —1 and is 
0{y~^~^‘^) as y —)• oo with —r + 2 < —1. In a similar fashion, we can also obtain 


c+oo 

Jwo 


dw 


[w^ taiih"‘{awzQ)]{l — F{w))dw > CwC(^ '^m{T) 


(3.113) 


with 


/ +00 

[2ytanh^(y) + 2y^tanh(y)(l — tanh^(y))] y~^'^~^Py, (3.114) 


for /3 sufficiently close to jdc- At this stage e > 0 is an arbitrary fixed quantity that will be chosen later 
(but independently of /3). By (j3.96p and (I3.97p . 


G{(3,z*Pz*,-^ - (1 - a^E[lF^]) ^ ^ G{(3,z*Pzr^ - (1 - a^E[lF^]) 

— ^0 — 


CwoF ^K{t) 

which, substituted in (|3.111l) and (13.1131) . gives 

„2TnUw2n"^-(^) / „2 


GwcF ^/c(r) 


(3.115) 


G{ld,ZQ)zl “ (1 “ ^ taxFp{awZ q)]{ 1 - F{w))dw 


<G{P,zl)zl 




k{T) 


-(l-a^E[lF^])- 


k{T) 


(3.116) 


22 






























By definition of D2{I3), 


G{/3,z*,)z*, 


K{t) 


+{l-a^K[W‘^]) 1 - 


m{T) 

'W) 


<D2i/3)<Gi^,z*o)z*o 


-i Mjr) 

k{T) 


+{l-a^W.[W‘^]) 1 - 


M(r) 

k{T) 


In the last step of the proof, we show that > 1- This can be done by properly choosing the arbitrary 
quantity e in (I3.114D . We will prove the first inequality, the second one can be obtained in the same way. 
Starting from (13.951) and (I3.114p . we introduce the functions Kb^r) and ma{T) for a > 0, 6 > 0 as 


/ + 00 j p+oo J 

— [y^ - ytanh(y)] ma(r) := J — [y^tanh^Py)] dy, (3.117) 

which coincide with K[t) and m{T) for 6 = 0 and a = e, respectively. By applying the integration by 
parts formula the two functions can be written as 


^+oo 

Kb{r) = -6^"^ + 6^"^tanh(6) + (r - 1) / (y^ - ytanh(y))y 

Jb 


"dy, 


ma(r) = —atanh(a) + (r — 1 


r+oo 

V ^ 


y‘^tanh^{y)y '^dy. 


Since 


mair) moir) ff^°°y‘^tanh.^{y)y '^dy 

iim , = -rr-^ -^- > 1, 


6^0+ 


0 + Kbir) K{t) [y2 - ytanh(y)] y-^dy 


(3.118) 

(3.119) 

(3.120) 


where the inequality can be proved by observing that y^ tanh^(y) > y^ — y tanh(y) for all y > 0, then for 
any e > 0 sufficiently small, 

m{T) 


K{t) 


> 1 . 


(3.121) 


Since G(/3,2 ;q)zq = 0{zq) = 0((/3 — f^cY) with s = > 1 and (1 — a^E[IT^]) = 0(/3 — /3c), with 

1 — a^E[IT^] < 0 for /3 > /3c and close to /3c (see (I3.72p i. we conclude that 0 < D2{I3) = 0{I3 — fdc), for the 
same values of /3. This proves that 

0<^ = O((/3-/3c)-i). (3.122) 

The previous equation together with zY^~^ = 0(/3 — /3c) allows us to conclude the proof along the same 
lines of the case with EpW^] < oo. Indeed, the upper bound (I3.77P is still valid in the present case, since 
only the first moment of W is involved. For the lower bound we argue as follows. Since, 1 — tanh^(a:) > 
1 — tanh(x) > 1 — X for x > 0, we have 


xi/3,0) > E 


r)?* 


rdzX 


.dzX 


= 1 - v^sinh(/3)E[IT]zo + Vsinh(/3)E[IT]-^ - smh{l3)vzl^ 

= 1 - Vsinh(/3)E[IT]0((/3 - + A/sinh(/3)E[IT]0((^ - ^c)"^) 

- sinh(/3)^.0((/3 - /3 c)'/("- 3))0((/3 - /3c)-'). 

The inequalities (j3.77p and (I3.123P imply (I3.79P concluding the proof of the theorem. 


(3.123) 

□ 
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4 Non-classical limit theorems at criticality: proof of Theorem 12.15 


In this section we prove Theorem 12.151 For this, we follow the strategy of the proof for the Curie-Weiss 
model (see e.g. m Theorem V.9.5]). It suffices to prove that for any real number r 


r, f ( \\ J^^expirz- fiz))dz 

(-» = /!L exp (-/(.)) d. ■ 


(4,1) 


As observed in |18| , the measure is approximately equal to the inhomogeneous Curie-Weiss measure 


1 




Wj^Wj 

ijG[N] 


1 




1 sinh /3 

cr)e^ ^iv 


2 (SliglAr] li’iO'i) 


,5®inh/3Ei, 

^ (TGOiv crGOiv 

where g{(j) is any bounded function defined in Ojv and Zjv is the associated normalization factor, i.e., 

1 sinh 8 ( v-' 

(Ei6[Ar] 


(4.2) 


Zj, = e" "Jv 


(4.3) 


We first prove the theorem for this measure Pjv, which is the rank-1 inhomogeneous Curie-Weiss model 
with (5 replaced with sinh(/3). 

For this, we use the Hubbard-Stratonovich identity to rewrite Pjv^exp as a fraction of two 

integrals of an exponential function in Lemma 14.11 in Section 14.11 We next split the analysis into the cases 
E[VF^] < oo and r E (3, 5) in Sections 14.21 and 14.31 respectively. For both these cases we analyze the 
exponents in the integrals and use Taylor expansions to show that they converge in Lemmas 14.21 and 14.41 
respectively. We then use dominated convergence to show that the integrals also converge in Lemmas 14.31 
and 14.51 respectively. The tail behavior of f{x) for r E (3,5) is analyzed in Lemma [4.61 Combining these 
results we conclude the proof of Theorem 12.151 in Section 14.41 we first prove the theorem for Pjv and then 
we show that the theorem also holds for Pjv in Lemma 14.71 Finally, in Section 14.51 we discuss how to 
adapt the proof to obtain the results on the scaling window. 


4.1 Rewrite of the moment generating function 


To ease the notation we first rescale S'jv by and later set A = 5/{8 + 1). We rewrite Pjv 
in the following lemma: 




Lemma 4.1 (Moment generating function of Sj^/N^). For B = 0, 


where 


with 


Pn 


exp 



J —oo 

e-^G'jvC;0)d^’ 

J—OO 


Gjv(z;r} = -z^ -E 


log cosh (q:jv(/3)IF)v5; + 


r \ 


N^/1 ’ 




/ 


sinh /? 

E[IF^] ■ 


(4.4) 


(4.5) 


(4.6) 
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Proof. We use the Hubbard-Stratonovich identity, i.e., we write e*^/^ = IE [e*^], with Z standard Gaussian, 
to obtain 


ZnPn ( exp (r 


r 1 sinh/S / 

Z^ie[JV] \ZielN] 


■ ZIie[JV] LV {Y.i!=[N] 


= 2^e[ n cosh(J^u;,Z + ^)] 

ielAT] ^ ^ 

We rewrite the sum in the exponential, using the fact that Wjv = where ?7jv is a uniformly chosen 

vertex in [N], as 


WiZ + ^] =2'"E e 


^Ei6[iv] log cosh 


ZpfPff exp r 


’-^)) = 2^E expjA^E log cosh 


I sinh (5 

NE[Wj^ 


W^Z + 




•)N /•oo 


J-c 


^/^expjA^E log cosh ^ajv(/3)bFiv—^ +jdz. (4.8) 


By substituting zj^pN for z, we get 


Zj.,Pj., ( exp (r 


= 2^Y^^^ J e ^^^/^expjA^E log cosh ^ajv(/3)W 


'iv^ + T7T rdz 


= r 


In a similar way we can rewrite 




(4.10) 


so that the lemma follows. 

4.2 Convergence for E[1W^] < oo 
We analyze the asymptotics of the function Gpf{z]r): 

Lemma 4.2 (Asymptotics of for E[W^] < oo). For j3 = j3c,N, B = 0 and E[W^] < oo, 

,r)- zry ^ ^uElW^f 
Proof. Taylor expanding logcosh(x) about x = 0 gives that 


(4.11) 


logcosh(x) = ^+ 0(x®). 


(4.12) 


We want to use this to analyze NGrf{z/N^P-,r) and hence need to analyze the second, fourth and sixth 


moment of 
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The second moment equals, using that A = S/(S + 1) = 3/4, 


E 


Z T \ ^~\ Z / ZT T 

J = sinh/3e,jvJ^jv;^ + 2^smh(^ cnHW^]— + 


zr E[Wj,] 

y/N NV 


+ o{l/N), 


(4.13) 


where we have used that sinh/3c,Ar = 1/i^jv in the second equality. 

For the fourth moment we use that by assumption the first four moments of are 0{1). Hence, for 
all r. 


E 


CiAr(/3c,7v)bFjv 


+ 


-V 


iVi/4 ' N^J \ E[lFjv] 

^ E[W^] z^ 

E[1T2]2 ]y 


!!hh^Eiir‘i-+ of— 

EIM/J2 ^ Wv-i 


+ 


1 


+ 


1 


+ 


1 


^^3/4+A ^^2/4+2A ^^l/4+3A JY'4A^ 

+ o(l/iV). (4.14) 


For the sixth moment, we have to be a bit more careful since E[VF®] is potentially infinite. We can, 
however, use that 


N N 

IE[VF®] = = {maxwi)^E[W, 

i=l i=l 


4] 


(4.15) 


i=l 
V 


It can easily be seen that max(^^ Wi = o{N^/^) when —> W and E[H/^] E[VF^] < oo. Hence, 


E 


Z \6- 


ajv(/3c,Af)bFjv^^y4 j 


sinh3/3c,7VTn,rw6i _ o(iVl/2)E[p^4] ^6 


E[W^ 


-E[w: 


In a similar way, it can be shown that 


E 


on{I3c,n)W^ 


' A^6/4 E[IF'2]3 iV6/4 

r \ " 


= o(l/iV). 


+ 


= o(l/iV). 


^Ni/i ' ata; j 

Putting everything together and using that the first four moments of Wj^ converge by assumption. 


(4.16) 


(4.17) 


lim NGi^{z/N^^*-,r) = lim [ — A^E 

N—^oo N—^oo \ 2 


( Z T \ 

aiv(/3c,iv)Wjv-^^ + 


/E[W] 1 E[W^] 4 

^12E[W2]2^- 


(4.18) 

□ 


From Lemma 14.21 it also follows that the integral converges as we show next: 

Lemma 4.3 (Convergence of the integral for E[IF^] < oo). For /3 = Pc,N, B = 0 and E[IF^] < oo, 


/ oo 1 ,4 

-oo J — i 


N^oo 


(4.19) 


' —CX> 
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Proof. We prove this lemma using dominated convergence. Hence, we need to find a lower bound on 
NGpf{z/N^^^;r). We first rewrite this function by using that 


E 




1 


2 VAri/4 


(4.20) 


Hence, 


1 




^ar^iPcN^W'^ - log cosh 

1 


G^{z/N^/^-,r) = E 

= E ^[aNif3c,N)WNj^'j - log cosh 


(4.21) 


Since 


-E 


d2 H 


log cosh {aM{l3c,N)WNj^ + - log cosh (^ajv(/3c,iv)lHiv-^^^ 


—— logcoshx) = 1 — (1 — tanm(x)) = tanm(x) > 0, 
dx^ 2 


the function ^x^ — log cosh x is convex and we can use Jensen’s inequality to bound 


E 


z \ 


aN{f3c,N)Wf,^^j - log cosh {aM{Pc,N)WNj^ 


> 2 [»NiPc,N)E[W^]j^j - log cosh {^a^{^c,N)E[WM]j^ 


Ie[Wj, 


iVl/4 


— log cosh 


/e[W^ 


AH/4 


As observed in the proof of m Theorem V.9.5], there exist positive constants A and e so that 


1 


x^ — log cosh X > d{x) := 


ex 

ex^ 


for |x| < A, 
for |x| > A. 


To bound the second term in (14.211) . we can use the Taylor expansion 

logcosh(a + x) = logcosh(a) + tanh(^)x, 
for some ^ G (a,a + x), and that | tanh(^)| < |^| < |a| + |x| to obtain 


E 


\ogcosh(^aNil3c,N)WNj^ + - log cosh (^ajv(/3c,Ar)Wj 




< E 


< E 


( Z T \ / Z 

aiv(/3c,iv)bHiv-^^^ + - log cosh (^ajv(/3c,Ar)bHiv-^^^ 

I I I I II 9 

Z ^ \ ” Z^ ^ 

aiv(/3c,iv)bHiv^^ = (^n{I3c,n)E[Wn] 


lE[Wr,] \zr\ 


+ 


N mG' 


Hence, 


,,-AfGjv(2/Afi/4.r) 


< exp 




^^^Azr\ + ^-Nd 


Ie[Wj, 


Ari/4 


)}. 


(4.22) 


(4.23) 


(4.24) 


(4.25) 


(4.26) 


(4.27) 
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which we use as the dominating function. Hence, we need to prove that the integral of this function over 
z G M is uniformly bounded. We split the integral as 


/I 




iVi/2 


Ie[Wm] z 


dz 




N 


1/4 


<A 


exp 




- \zr\ + —TTTT — JAd 


Ari/2 


'IE[W^ 




)} 


dz 


+ 


E[W^r] 


“^JV jvi/4 




exp 


{/ 


IE[Wiv], , , r2 


/E[W^] z 

I^JV iVl/4 


)} 


dz 


(4.28) 


The first integral equals 

/, 


E[Wjv] ^ 
I'AT Arl/4 


<A 


exp 




— -\zr\ + —rw: — £ 


E[W^]4 4 


A^V2 E[W2]2' 


)}d=, 


(4.29) 


which clearly is uniformly bounded. The second integral equals 

..2 


/ 


E[W^r] ,, 

vn ivl /4 


>A 


exp 


r /e[w^] 




(4.30) 


iVl/4 


/, 


E[W,v] „ 
^'N Arl/2 


>A 


exp 


{/ 


E[W^] |j/r| r2 E[W^j 2 




where we have substituted y = zA^^/*^. This converges to zero for N -A 00 , because the integral is uniformly 
bounded. 

Together with the pointwise convergence proved in Lemma 14.21 this proves Lemma [4.31 □ 

4.3 Convergence for re (3,5) 

We next analyze Gi^{z',r) for r E (3,5), assuming Condition 12.51 

Lemma 4.4 (Asymptotics of Gjv for r E (3,5)). Assume that Condition \2.5\h ii) holds. For (3 = fdc^N, 
B = 0 and r E (3, 5), 


lim NG^{z/N^^^^-^^;r) = -zrJ^^ + f (J 
TV —>00 y V \ V 


E[W] 


(4.31) 


where /(z) is defined in (I2.37p . 
Proof. Dehne the function 


j_ 2 

g{w, z) = - {aN{fic,N)wz) - log cosh {aN{fic,N)wz ), 


so that we can rewrite, in a similar way as in (14.211) . 
NGn (z/N^;r) = iVE[ 5 (Wjv,z/iV^)] 


(4.32) 

(4.33) 


- iVE 


log cosh 




— log cosh 


-.W, 
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By the definition of we can rewrite 


N 




(4.34) 


2 = 1 


With the deterministic choice of the weights as in (I2.8j) . 


1 


- log cosh 


C,nZ 


2 \^y^E[W2] 

From this it clearly follows that, for all f > 1, 


— log cosh 


An/l-r-l), 

1 CujZ 

TT 




(4.35) 


C,nZ 


lim g{wi, ^^) = “ ( — , —rn —TT 

N^oo ^ ' ^2 I y^E[W2] iV(^-l) 


— log cosh 


C,„z 


y^E[PF2] zV(r-l) 


(4.36) 


It remains to show that also the sum converges, which we do using dominated convergence. For this, we 
use a Taylor expansion of log cosh(x) about x = 0 up to the fourth order 


4 2 4 

X X X 


logcosh(x) = Y + (-2 + 2tanh^(^) + 6tanh^(,f)(l - tanh^(,f))) Y - Y ~ 12’ 
for some ^ E (0, x). Hence, 


g{wi,z/N^/<^--^^) < 


1 


1 


CwZ 


1 


(CwzY 


12 y y/E[iy2] fVH-i) y 12 E[iy2]2 + o(l) • 

Since r E (3, 5), it holds that 4/(r — 1) > 1, so that 

1 (cwzY 


(4.37) 


(4.38) 


7V^ ^ 12 E[W2]2 + o(l) H/H-1) 


< oo. 


(4.39) 


We conclude that 
N 

N—>oo 


lim '^g{wi,z/N^^^'^ = X] 

>/_i.rv-\ * ^ ^ 


C,nZ 


2=1 


2 = 1 


v/E[W2] H/H-i) 


— log cosh 


CwZ 


{-!- 

\v^E[W2] H/H-i) 


= / 


E[IF] 


(4.40) 


where in the last equality we have used that E[IF] = This is in turn a consequence of the following 

explicit computation giving an upper and lower bound on E[IF)v] matching in the limit N —)• oo. An upper 
bound on the first moment is given by 


1 ^ 


2=1 
r — 1 


= c„JV-Sl < c.JV-Sf (^1 + 1 " 


— C-’ii 


1 T-2 

- c,„ - -N ^-1, 


r-2 r-2 


(4.41) 
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and a lower bound by 




pN 

Ji 


From this it indeed follows that 


i V(^ — \-c^- ——1. 

r — z r — 2 

r — 1 


E[i^] = lim E[iF^] = -- 

N^oo T — 2 


To analyze the second term in (14.33^ . we can use the Taylor expansions 

log cosh(a + x) = log cosh(a) + tanh(a)x + (1 — tanh^(^))x^ 

= logcosh(a) + (a — tanhC(l — tanh^ C)a‘^)x + (1 — tanh^(^))x^, 


for some ^ G (a, a + x) and C € (0, a). This gives 
log cosh ( 


IVE 




V^EiM^ JV^ iVV 
= N 


log cosh 






Ie[Wj, 




- IVE 


+ iVE 


tanhC(l — tanm C)'^i 
^2 


r 




(1 — tanh^(^)) 


J iV2A 


Ie[w^ 


-zr + o(l), 


where the last equality follows from A = and r G (3, 5). 

Again it follows that also the integral converges: 

Lemma 4.5 (Convergence of the integral for r G (3,5)). For /3 = f3c,N, B = 0 and r G (3, 5), 


lim 

N^OO 




J —CO 


dz. 


Proof. We again start from the rewrite of in (14.211) . As before, 


Af 1 . 


C,nZ 




— log cosh 


C,,,z 


VEpTfl 


where it is easy to see that the summands are positive and decreasing in i. Hence 

pN .-1 , 1 ^ .V ^ 2 


iVE[i,(H/„, VW--.)] > - logcosh ( 

We want to use (14.241) , and hence split the integral in the region where | 


1 


c,„z 


(4.42) 

(4.43) 


(4.44) 


(4.45) 

□ 


(4.46) 


(4.47) 


dy. (4.48) 


V^EiW|]yV(-i). 

^nw%] I smaller 
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than A. This gives 


NK[g{W^,z/N—)]>e 


iE[Ty2] 


/ 


cw\z\ 




„4 ^4 pN 


dy 


= e 


4^2 r-1 / fA^E^Ur-3 


^\Wn] t - 3 I V Cu,\z\ 


- 1 


— £; 


C,„\z\ 


/ 5^ _ / ^v/Ep^ ^^-(5-r) 

■E[Ty2]5_T|^ V 

= yti|zr(^+i) - k2z‘^ - o{l)z^ + hlzl^-^, 
for the proper constants ki,k 2 ,k 3 > 0. Since 9 — r > 4, 

/ CO 

^-ki\z\-^^+^Ak2Z^+o(l)z^-k3\z\^-'^ 

-CO 


< OC. 


(4.49) 


(4.50) 


Together with the pointwise convergence in the previous lemma, this proves this lemma for r = 0. For 
r 7 ^ 0, the proof can be adapted as for the case E[Ty^] < oo. □ 

We next analyze the large x behavior of f{x) arising in the density of the limiting random variable: 

Lemma 4.6 (Asymptotics of / for r G (3, 5)). For r G (3,5), 


lim 

x—>-oo X 


fix) 


T —1 


{vzi) L (^ 


Proof. We first prove that the integral is hnite. For this, dehne 


Hy) = - log coshy. 


so that h{y) > 0. Then, 

f( 


— log cosh 


2y2/(r-l) 


l/(r-l) 


y 


f 

Jo 


dy= h 


1 


l/(r-l) 


dy. 


Since log coshy > 0, we have h{y) < ^y^, and hence 


h 


< 


y^/ir-l) ) - 2y2/(^-l) 


(4.52) 


(4.53) 


(4.54) 


This is integrable for y —>■ 0, because 2/(r — 1) < 1 for r G (3,5). 
Using (I4.37p . for y large, 


h 


1 


< 


1 1 


y^lP 1) j 12y^/U 1) 

This is integrable for y —>■ oo, because 4/(t — 1) > 1 for r G (3, 5). 

To prove that f(x)jx'^~^ converges to the integral as x —>■ oo we rewrite, with a = {r — 2)/(r — 1) 


(4.55) 


fix 


1 1 °° 


2 = 1 




= a 


T —1 


1 


{axY 


tE'* 

2=1 


1 


if {ax 


.T— 1 




= a 




(4.56) 

□ 
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4.4 Proof of Theorem 12.151 

We can now prove Theorem 12.151 for the measure Pm'- 


Proof of Theorem \2.15\ for the measure P^- We can do a change of variables so that 


f 


/•LX.) 


z. 


Hence, using Lemma HTT 


Pm ( exp (r 


5. 


TOO g-ArG^(z/iVi/(«+i);r.)^^ 

J —oo 


]^6/{6+l))) TOO g-7VGjvWAfi/i‘^+i);0)(i_2 ’ 

J— OO 

It follows from Lemma 14.31 for E[IT^] < oo and from Lemma [4.51 for r G (3,5) that 


lim Pm ( exp ( r 


N^OO 


N^/{S+i) 


roo 

J-oo^ 


dz f°° e^^-^(^)da 

J —CO 


roo 

J-oo^ 




dz 


f°° e ^(*)dx ’ 

J —OO 


(4.57) 


(4.58) 


(4.59) 


where we made the change of variables x = \ in both integrals to obtain the last equality. 

As mentioned, this is sufficient to prove the convergence in distribution of jygf(a+i) to the random 
variable X (see [HI Theorem A.8.7(a)]). 

For the case E[VF^] < oo, 


1 E[wd 4 

lim M. = lim 1 ] 


x^oo X 


1+5 


12E[IT]^ 


For T G (3, 5), the proof that lim 3 ;^oo = C is given in Lemma IT 


(4.60) 


□ 


It remains to show that the statement of Theorem 12.151 also holds for the measure Pm- This follows 
from the following lemma: 

Lemma 4.7. For E[VF^] < oo and r G (3, 5), 


hlf)) - “p hdt j 1=“■ 




(4.61) 


Proof. As shown in [T8 


where 


Jij = - log 


PN{g{(T)) = 


1, f e^pij + {I- pij) 




9(o')e 


-Y ■ 
2 


2 


E 




g2 Si,i6[A''] 


e ^Pij + {l-Pij) 


= pij sinh fd — pfj sinh /3(cosh /3 — 1) + 0{pfj), 


(4.62) 


(4.63) 


where we have used the Taylor expansion of log(l + x) about x = 0 in the last equality. Hence, using (12.5L 


g2 ^i,2e[JV] ' 


JijaiCTj ^ I Ei,j6[iV] [[ eJ+ZUj -^) sinh/l(cosh/l-l)+0(p?.)y,aj lsinh/3Eij 


e[^l f-M 


1 Wj^Wj 

= ■ QEM{cr)^2^^^^hT,iJelN] ^ 


(4.64) 
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Hence, we can rewrite (14.621) as 


PN{g{(T)) = 


/ ^(tG^/v 




(4.65) 


Combining this with the Cauchy-Schwarz inequality gives 




;’A.(exp (r^) - P»(e'^»<">))) 


< 


Pjv ( exp ( 2r 


P^(eEivW) 

P^(e®ivW) 


Pjv ( exp ( 2r 


. SiV 


p^(e2EivW _p^(ei?iv(<T))^ 




(4.66) 


From (14.591) . it follows that the first square root converges as N ^ oo. We next analyze Ej^[a) and show 
that E{^{a) —>• 0 in probability w.r.t. Pjv- We also show that Ej^[a) is uniformly bounded from above, so 
that the lemma follows by dominated convergence. 

We first analyze the contribution of the 0(p? ) terms in Ef^{a). Note that 


E Pij = Y1 

i,j&[N] i,j&[N] 


WiWj 


+ WiWj 


< 


E 

i,je[N] 


WiWj 

~J7 


E 


£3 

^ Ee[Ar] 


Wi 


(4.67) 


For E[VF^] —)• E[VF^] < oo it holds that max* Wi = o{'/N). Hence, 


'ielAf] 

because = 0{N). Hence 

= \ Y1 

i,j£[N] 

= -- E 

2 ^ 

i,je[N] 


wf^ < (maxu;*)^^ Y 


(4.68) 


WiWj WiWj 


Irq+WiWj if, 


wfw'j 

^Jv(^JV + WiWj 


'ieiw] 


sinh (5 — pfj sinh /3(cosh /3 — 1) + 0{pfj) ) a^aj 


■ sinh fS + 


WiWj 


£n + WiWj 


sinh /3(cosh /3 — 1) I crjCTj + o(l) 


--sinh/3 cosh/3 ^ —j^cTjCTj + o(l) 


wfw‘j 


£2 

i,je[N] ^ 


/ 2 \ 2 
sinh/?cosh/3[ E 


(4.69) 


where the third equality can be proved as in the analysis ofp?-. Hence, Ej^{a) is indeed uniformly bounded 
from above, so that is uniformly bounded. 
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It remains to prove that —)• 0 in probability w.r.t. Pjv- We define IV = ^<^0 so that 


Epf = — - sinh/3cosh/3 + o(l). 


(4.70) 


We analyze the moment generating function of Yjv the same way as S^/N^. That is, we use the Hubbard- 
Stratonovich identity to rewrite 

1 sinh /3 _ \2 

V prYNp2^^(E,ie[N]^i'^i) 

^iv fe ) — - 


E 




1 Sinn p /T—' \9 

g2-5]/Ci6[iV] 


(4.71) 


E 


ctGO.n 


E 




g _.e[JV] JVE[W^]^* + V iVE[Wjv] ^i6[JV] 


sinh (3 




E 


/ sinh /3 _ ^7 

Jg p\l JVE/Wjvl ^ie[JV] 

^iVE[logcosh I z; 


E 


^7Ve[ log cosh (y^^M/^z) I Z] 


(4.72) 


r°° g 2 :^/ 2 +A''e[ log cosh (^jvEilErl~*~\/^l'V^ltt"jv 2 :)] 

_ J — CO _ 

coo g-P/2+AfE[ log cosh 

J— CO 

We do a change of variables replacing zjy/N by z, so that 

_ r°° Af5;^/2+AfE[logcosh jvElWjyl ~*~'\/ElWj^l 

Pn ^_ 

J —CO 

CCX. ^-7VG^(^;0)+7Ve[ log cosh + -log cosh (y|^14^^.)] 

_ J— CO 

“ f°° e-^G^(^;0)dz 

J—CO 

coo ^-7VG^(V7VV(^+i);0)+7VE[logcosh + -log cosh 

_ J —CO 

“ r°° g-iVGjv(2/Afi/(‘^+i);0)(i2 ’ 

J —CO 

where we did another change of variable in the last equality. 

In Lemmas 14.21 and 14.4[ we proved that NG!^(z/N^^^^~^^^;0) converges for /3 = (3c- We Taylor expand 
the remaining term, 


A^E 


' sinh (3 . 




= E 




(4.73) 


ForE[W3] ^E[IT3 ] < 00 , which includes power-law distributions with r > 4, we can use that | tanh(x)| < 
|x|, so that 



r f / 

E 

tanh ( .. / 


sinh (3 






E[ITjv] ^iVi/(<5+i)y E[IT^]J - y E[IT^] A^V(^+i) E[W 


< 


sinh/3 \zr\ E[IT^] 


= 0 ( 1 ). 


(4.74) 
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For T G (3,4] we use the deterministic choice of the weights as in (12.Sp and d = r — 2 to rewrite 

, N 



r / / 

E 

tanh ( .. / 


sinh/3^^^ z 
Wn- 




IE[hFjv] iVi/(<5+i) J E[hFj 

1 


tanh 


sinh /3 z 

rWi- 


2=1 

„2 / aT\ 2/(r-l) 


iE[hFjv] *ivi/(r-i); 


Wt 


t , / /sinh/3 Cu,z \ rcl, fN\‘ 

N ^ VY E[VF^] iVF-i) y* E[fF^] V i J 

^2/(r-l)-l ^2/(r-l)-l / sinh/3 |rz|c3 ^ 3/(r-l) 


< Flci; 


E[hFj 


For T > 3 the first term is o(l). For r € (3,4), 


(4.75) 


^ pN 1 

^2/(r-l)-l ^ .-S/Cr-l) < J i-3/{r-l)^- ^ ^_Z_ _ AT-VF-l)^ ^ 


i=2 

whereas for r = 4 


^2/(r-l)-l ^ ^-3/(r-l) < / r3/(^-l)df = iV"^ log N = o(l). 

i=2 


(4.77) 


Hence, in all cases the integrands in the numerator and denominator of (14.7211 have the same limit. In 
Lemmas 14.31 and 14.51 it is proved that the integral in the denominator converges. Since 



r f / 

E 

tanh ( .. / 


sinh (3 z 

W N- 


W^. 


< 


rE[LF2 

E[LF^; 


= 0 ( 1 ), 


(4.78) 


E[LFjv] ^iVi/(<5+i); e[LFjv 
it follows by dominated convergence that the integral in the numerator has the same limit. Hence, 

lim Pjv (e’'^^) = 1, (4.79) 

N—>-oo ^ ' 


from which it follows that 0 in probability w.r.t. P^. Hence, also — ^ sinh/? cosh/3 —?> 0 in 

probability w.r.t. Pjv- Since o(l) also converges to 0 in probability, so does the sum: 

Fijv = sinh/3cosh/3 Y^^ + o(l) —;■ 0 in probability w.r.t. Pjv- (4.80) 

□ 


Remark 4.8 (Sharp asymptotics of the partition function). It follows from the changes of variables 
in (33]) QjThd (|4 . 57p iyhojt 

Z^iPc, 0) = Hiv02+i/(<5+i)2Af(i + o(l)). 

For E[VF^] < oo, this exponent equals 1/2 + l/((5 + 1) = 3/4, whereas for r G (3, 5), it is l/2 + l/{5 + 1) = 
(r + l)/(2r — 2). Thus the partition function has finite-size power-law corrections (in agreement with 
where the classical Curie-Weiss model is considered). 
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4.5 Scaling window 


Instead of looking at the inverse temperature sequence /3jv = /3c,iv we can also look at ~\~b/N 

for some constant b. The analysis still works and the limiting density instead becomes 


exp <J ^ cosh(/?c)®^x2 - f{x) } . 


(4.81) 


To see why this is correct we look at the following second moment, which shows up in the expansion of 
Gjv, see (I4.13P : 


-E 

2 


' sinh(/3c,Ar + b/N ^+i ^ 


■IT, 


r \ 


2i 




nw^] 

^-1 E[IT?rl / a -1 zr 

sinh(/3c,Ar + + Vsinh(/3c,Ar + b/N^)E[Wr,]- + oil/N). 


(4.82) 


2iV2/(5+l) 

In the first term, we Taylor expand the sine hyperbolic about /3c, at, which gives 
sinh(/3c,Ar + b/N'^) = sinh(/3c,7v) + cosh (/3c, at)6/+ O 
For the other term, and also for the other terms in the expansion of G^, it suffices to note that 


Hence, (|4.82p equals 


smr 


+ 


HPc,n + b/Ns+i) = Jsmh{/3c,N) + 0 (l//V<5+i 


(4.83) 


(4.84) 


2iV2/(^+l) 2A^2/((5+1) 

-2 6^2 


cosh(/3c,Ar) 


b E[W^] zr K[Wj 


+ 


jYlfiE[IT^] ' N ,/e^- 
E[IT2] zr E[IT,v] 


+ o(l/iV) 


2iVV(^+i) + 2jv + iV Tlpf + 


(4.85) 


In the expansion of the first term in (14.8511 drops as usual, whereas the second term in 

(I4.85P remains. After multiplication by N (cf. (I4.18p b one has 




Va/eIif; 


0 ( 1 ) 


(4.86) 


Using the substitution x = the above converges in the limit A" —)• oo to the exponent in (I4.8ip . 

as required. □ 


Limit distribution at /3c instead of /3c, at* In the above, we look at the inverse temperature sequence 
Pn = /3c,Af and then take the limit N —>• oo. Alternatively, we could immediately start with j3 = /3c. The 
scaling limit that will be seen depends on the speed at which approaches v. Indeed, from (I2.24p and 
(I2.34p . one has /3c — (3c ,n = 0{i' — v^)- 

We investigate this for the deterministic weights according to (12.Sp . and first investigate how close 
j/jv is to Tz. By |3l Lemma 2.2], = v + C,N~'^ + o{N~^) with rj = (t — 3)/(t — 1) and C an explicit 

non-zero constant. Thus, for r > 5, = r' + Hence, the results stay the same (see the previous 

discussion). 

a-i i-i 

When r G (3,5), instead, = v + CN~‘^ + o{N~'^) = u + CN «+i -|-o(A <5+1), so we are shifted inside 
the critical window (see the previous discussion). Hence, in this case the limiting distribution changes. □ 
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Table 1: List of symbols used 


Symbol 

Definition 

Description 

N 


Number of vertices 

[N] 

,Ar} 

Set of first N positive integers 

/? 


Inverse temperature 

B 


External field 



Hamiltonian 

Zn 


Partition function 


limAT^oo jj log Zn 

Pressure of inhomogeneous Curie-Weiss model 

Wi 


Weight of vertex i 

w 

{wi,.. .,wn) 

Sequence of weights 

In 

\~^N 

l^i=i Wi 

Total weight 

GRGn{w) 


Generalized random graph with weights w and N vertices 

Pij 

WiWj 

Probability of an edge between vertices i and j in GRGn/w) 

iN+WiWj 

Wn 


Weight of uniformly chosen vertex 

w 


v> 

Random variable chosen such that Wn — W 

y 

K[W‘^]/¥.[W] 

Size-biased weight 

Vn 

]E[1T2]/E[1T^] 

Its finite volume analogue 

T 


Power-law exponent 

Qn 


Law of the random graphs 

Pn 


Annealed Ising measure 

'tpN 

log Qn{Zn) 

Annealed pressure 

Sn 


Total spin 

Mn 

Pn/Sn/N) 

Annealed magnetization 

Xn 

-^Mn 

Annealed susceptibility 



limAT^oo of i/jn^Mn^Xn, respectively 

/3c 

asinli(l/z^) 

Annealed critical inverse temperature 

13c, N 

asinli(l/z^jv) 

Its finite volume analogue 

/3,<5,7,7' 


Critical exponents, see Def. 12.91 

z* 


Pixed point of (]2.23|) 

Pn 


Curie-Weiss approximation of Pn, see (14.2p 

Zn 


Curie-Weiss approximation of Zn, see (14.3p 


37 









Acknowledgments. We thank Aernout van Enter for helpful discussions on inhomogeneous versions of the Curie- 
Weiss models. We thank Institute Henri Poincare for the hospitality during the trimester “Disordered systems, 
random spatial processes and their applications”. We acknowledge financial support from the Italian Research 
Funding Agency (MIUR) through EIRE project grant n. RBFR10N90W. The work of RvdH is supported in part 
by the Netherlands Organisation for Scientific Research (NWO) through VICI grant 639.033.806 and the Gravitation 
Networks grant 024.002.003. 

References 

[11 R. Albert and A.-L. Barabasi, Statistical mechanics of complex networks. Review of Modern Physics, 
74, 47-97 (2002). 

[2] I. Benjamini and O. Schramm, Recurrence of distributional limits of finite planar graphs. Selected Works 
of Oded Schramm, 533-545, Springer New York (2011). 

[3] S. Bhamidi, R. van der Hofstad, and J. S. H. van Leeuwaarden, Novel scaling limits for critical 
inhomogeneous random graphs. The Annals of Probability, 40 (6), 299-2361 (2012). 

[4] G. Bianconi, Superconductor-insulator transition on annealed complex networks. Physical Review E, 85(6), 
061113 (2012). 

[5] B. Bollobas, S. Janson, and O. Riordan, The phase transition in inhomogeneous random graphs. 
Random Structures & Algorithms, 31(1), 3-122 (2007). 

[6] A. Bovier, Statistical mechanics of disordered systems: a mathematical perspective. Cambridge University 
Press (2006). 

[7] F. Comets, N. Gantert and O. Zeitouni, Quenched, annealed and functional large deviations for one¬ 
dimensional random walk in random environment. Probability Theory and Related Fields, 118(1), 65-114 
( 2000 ). 

[8] A. Dembo and a. Montanari, Ising models on locally tree-like graphs. Annals of Applied Probability, 20, 
565-592 (2010). 

[9] A. Dembo and A. Montanari, Gibbs measures and phase transitions on sparse random graphs. Brazilian 
Journal of Probability and Statistics, 24, 137-211 (2010). 

[10] S. Dommers, C. GiardinA and R. van der Hofstad, Ising models on power-law random graphs. Journal 
of Statistical Physics, 141(4), 638-660 (2010). 

[11] S. Dommers, C. GiardinA and R. van der Hofstad, Ising critical exponents on random trees and 
graphs. Communications in Mathematical Physics, 328(1), 355-395 (2014). 

[12] S.N. Dorogovtsev, A.V. Goltsev and J.F.F. Mendes, Ising models on networks with an arbitrary 
distribution of connections. Physical Review E, 66, 016104 (2002). 

[13] S.N. Dorogovtsev, A.V. Goltsev and J.F.F. Mendes, Critical phenomena in complex networks. 
Reviews of Modern Physics, 80(4), 1275-1335 (2008). 

[14] R.S. Ellis, Entropy, Large Deviations, and Statistical Mechanics. Springer-Verlag, New York (1985). 

[15] R.S. Ellis, C.M. Newman, Limit theorems for sums of dependent random variables occurring in statistical 
mechanics. Zeitschrift fiir Wahrscheinlichkeitstheorie und verwandte Gebiete, 44(2), 117-139 (1978). 

[16] R.S. Ellis, C.M. Newman, The statistics of Curie-Weiss models. Journal of Statistical Physics, 19(2), 
149-161 (1978). 


38 



[17] C. Giardina, C. Giberti, R. van der Hofstad and M.L. Prioriello, Quenched central limit theorems 
for the Ising model on random graphs. Journal of Statistical Physics, 160, 1623-1657 (2015). 

[18] G. Giardina, G. Giberti, R. van der Hofstad and M.L. Prioriello, Annealed central limit theorems 
for the Ising model on random graphs. ALEA, Latin American Journal of Probability and Mathematical 
Statistics, 13(1), 121-161 (2016). 

[19] R. VAN DER Hofstad, Random graphs and complex networks. Vol. I. Lecture notes. Preprint, (2014). 

]20] M. Krasnytska, B. Berche, Y. Holovatch and R. Kenna, Violation of Lee-Yang circle theorem for 
Ising phase transitions on complex networks. Europhysics Letters, 111(6), 60009 (2015). 

[21] M. Leone, A. Vazquez, A. Vespignani and R. Zecchina, Ferromagnetic ordering in graphs with 
arbitrary degree distribution. The European Physical Journal B, 28, 191-197 (2002). 

[22] A. Montanari, E. Mossel and A. Sly, The weak limit of Ising models on locally tree-like graphs. 
Probability Theory and Related Fields, 152, 31-51 (2012). 


39 



